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ABSTRACT 

We study the effects of halo substructure and a distribution in the concentration parameter of 
haloes on large-scale structure statistics. The effects on the power spectrum and bispectrum 
are studied on the smallest scales accessible from future surveys. We compare halo-model 
predictions with results based on A^-body simulations, but also extend our predictions to 10- 
kpc scales which will be probed by future simulations. We find that weak-lensing surveys 
proposed for the coming decade can probe the power spectrum on small enough scales to de- 
tect substructure in massive haloes. We discuss the prospects of constraining the mass fraction 
in substructure in view of partial degeneracies with parameters such as the tilt and running of 
the primordial power spectrum. 

Key words: galaxies: haloes - cosmological parameters - dark matter. 



1 INTRODUCTION 

The substructure of dark-matter haloes is a topic of interest in 
cosmology, especially due to its relevance to the cold-dark-matter 
model, which requires some degree of dumpiness to the matter dis- 
tribution. In a recent paper, Kochanek & Dalai (2003) argue that 
substructure provides the best explanation (originally proposed by 
Mao & Schneider 1998) of the anomalous flux ratios found in grav- 
itational lens systems. However, the observed number of dwarf 
galaxy satellites in the Local Group is more than an order of mag- 
nitude smal ler than the number of subhaloes foun d in numerical 
simulations < KlvDin et alll999l:lMoore et alll9993) . 

Observable statistics of the matter distribution provide an ap- 
proach to substructure studies that is complementary to studies of 
strong-lens systems or direct observations of satellite galaxies. Sub- 
structure effects are expected to produce an observable effect in the 
power spectrum in the next generation of lensing surveys. As we 
show in this paper, substructure effects in the weak-lensing regime 
probe cluster scales primarily, as opposed to strong-lensing sys- 
tems, which probe galactic scales. 

The halo model provides a description of non-linear grav- 
itational clustering which allows statistical quantities of interest 
to be efficiently co mputed. The essential c oncepts of the halo 
model were used by iNevman & ScottI il952h to describe the spa- 
tial distribution of galaxies over fifty ye ars ago. The ir results were 
later generalized by Scherrer & Bertschi nged il99lh . By combin- 
ing the formalism of Scherrer & Bertschinger with results from 
A^-body simulations such as a halo density profile (e.g. that of 
iNavarro. Frenk & White||l 997h and a hal o mass function and bias 
parameters (e.g. Ishetl^fe TormenI ^9^, one is able to calcu- 
late various dark-matter and galaxy clustering statistics <Ma & Fr\< 
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prehensive review of the halo model. 

Most work on the halo model approximates haloes as spher- 
ically symmetric objects, described by a mass density which is a 
smoothly-decreasing function only of the distance from the halo 
centre. However, the density distribution of haloes which form 
in A^-body simulat ions is not at all smooth iNavarro et aLT997l : 
iMoore et alll999a) . Rather, one fi nds that about 10 per cent of the 
mass is associated with subhaloes iTormen. Diaferio & SveJl998t 
Ghigna et al. 2000). In the large-scale structure literature to date, 
the effects of halo substructure are typically assumed to be small 
for purposes of calculating statistics of the dark-matter distribu- 
tion. This paper explores the validity of this assumption. We will 
make quantitative statements about the size of substructure effects 
on the 2D and 3D dark-matter power spectra and bispectra. We will 
identify the length and angular scales at which halo substructure 
becomes important. 

The results of A'^-body simul ations also suggest a distribution 
of halo concentration parameters ( Bullock et al. 2001; Jing 2000). 
Like substructure, this effect is often ignored, thou gh the formal- 
ism to include it in the halo model is straightforward iCoorav & Hii 
I2OOII) . This paper will quantify the effect of a concentration param- 
eter distribution on the power spectra. 



2 THEORY 

To model substructure in a halo, the total mass A/ of the halo is di- 
vided into a smoothly distributed component of mass A/s, and the 
remaining mass m is placed in subhaloes. We will use the terms 
subhalo and subclump interchangeably. Ignoring substructure, the 
halo model requires specification of the halo density profile and the 
mass function, i.e. the number density of haloes of a given mass. 
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To include substructure, one must specify not only the density 
profile of the smooth component, the Fourier transform of which 
we denote Ms U(k\ Ms), but also the density profile of the sub- 
haloes, with Fourier transform mu{k; m). The mass factors sim- 
ply provide the proper normalization; U{k; Ms) and u{k; m) are 
normalized to unity. Suitable mass functions for the number den- 
sity both of parent haloes, dN/dM, and of subhaloes, dn/dm, 
are also needed. In addition, the distribution of subhaloes within 
the parent halo must be specified. We denote its Fourier tr ansform 
by Uc{k; Ms). Using these ineredients. lSheth & JainI i2003 ^ derive 
expressions for the power spectrum and the bispectrum. We give 
some results from their paper that are relevant to this work in the 
Appendix. 

We use the mass function of ISheth & Tormer] ^99^ for the 
number density of parent haloes of a given mass M: 



dM 



dN{M) 



■A[l + {av)~'']^e 



-av/2 
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dM M ^ - V 

where p denotes the average density of the Universe at the present 
epoch, and 

1 2 



v{M, z) ■- 



D{z)a{M) 



(2) 



is the peak height defined in terms of 5c{z), the threshold overden- 
sity in the spherical-collapse model, D{z), the growth factor for 
density perturbations, and a{AI), the rms fluctuations in the matter 
density field at the present epoch smoothed with a top-hat window 
of radius (3M/47tp)^''^. The normalization constant is fixed by re- 
quiring that all mass be contained within a parent halo, an d turns 
out to be ^ = 0.116. For Sc, we use the fitting formula of lHenrvl 
i200(]|) . We use haloes with mean interior density 180 times the 
background density and ta ke a = 0.67 and p = 0.33, so that our 
mass function is universal <Whitel2002h . The notation of this paper 
reflects our choice of halo boundary: r := risob, M :— Misob and 
c := risob/fs, where rs is the usual scale radius parameter of the 
halo density profile. 

We use the density profile of iNavarro et alj il997b . hereafter 
NFW, for both the smooth component and the subhaloes. It's 
Fourier transform is most efficient to compute using sine and co- 
sine integrals: 



U{k, m) = cos K {Ci + c)] - Ci(K)} 
with K := krisob/c, 

j-x 

Si{x) - I dr ■ 
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and 



Ci(x) 



dr ■ 
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We assume the distribution of subhaloes within a parent halo fol- 
lows the mass distribution of the smooth component. That is, 
Uc{k; Ms) is also given by the NFW profile, normalized to unity. 
Simulations suggest that the abundance of subhaloes in the par- 
ent halo traces the density profile of the parent quite well except 
in the most dense region near its centre where the abundance of 
subhaloes is lessened due to tidal disruptions. The subhalo abun- 
dance in the inner region is not yet well resolved in simulations. 
IChen. Kravtsov & KeetonI i2003') propose a modified NFW profile 
that is constant or zero inside of a core radius which they take to 
be 10-20 kpc for a lO^^-M© halo. For a halo of this mass, only 4 
per cent of its matter lies within the core radius. Hence, were we to 



use the latter profile, the effect would be to move only 4 per cent 
of our subhaloes further from the centre. We will comment briefly 
on the effect we expect this to have on the power spectrum in the 
Discussion. 

W e use the concentration parameters found bv lBullock et alJ 
i200ll) . for all three profiles at redshift z, which are 



c{M,z) = 9(1 + 



Mvir(M) 



M, 



for the smooth component of a halo and 



c(m, z) = 7.5 {1 + z) ^ 



mvii{m) 



M* 



(smooth) , 



(subhalo), 



(6) 



(7) 



for all subhaloes. The non-linear mass scale AI, is defined by the 
relation !/(Af, , 0) = 1. The translation Afvir(Af) b etween different 
halo mass definitions is derived in the appendix of lHu & Kravtsov) 
( 2003). Bullock et al. find that the concentration parameter for sub- 
haloes has some dependence on the local density around the sub- 
halo. We ignore this effect here, as it is difficult to include in our 
model, and we do not believe inclusion wo uld seriously affect our 
conclusions. 'H uffenbereer & Seliakli2003l) find a different concen- 
tration parame ter by minim ization against the non-linear fitting 
formula of lSmith et aljj2003h : 



c{M, z) = 11 (1 + z)-' 



(8) 



Numerical simulations suggest a log-normal distribution for 
the concentration parameter about the mean value c. 



p(c; m, z) dc - 



d(lnc) 



: exp 



In^ [c/c(m, z)] 



2<. 



(9) 



with a width ai„^ = 0.3 <Jindl2000HBullock et alj|200ll) . Unless 
otherwise noted, we use the mean values given in equations ^6) and 

o 

We assume a biasing prescription for the clustering of parent 
haloes. That is, the halo power spectrum is given in terms of the 
linear mass power spectrum Pl (k) as 



P{k,z;Mi,M2) = bi{ui)bi{iy2)PLik,z), 



(10) 



whe re Uj := v(M i,z). We use the bias parameter of 
Shet h & Tormed <1999h : 



av — 1 



2p/4 
1 + {av)p' 



For the subhaloes, we assume a power-law mass function su; 
geste d by numerical simulations ( Tor men et al...l998 : Ghigna et 



(11) 



dn{m\M) 
dm 
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ml m 



(12) 



where p < 1. The normalization A*'o is determined by the mass 
fraction contained in subhaloes, /: 



/ 



M -Ms 
M 



dm 



m dn(m; M) 
M dm 



(13) 



Numerical simulations suggest p « 0.9 and / ~ 10 per cent. We 
use these values unless otherwise noted. In addition, simulations 
suggest that subhaloes heavier than 1 per cent of the total mas s in 
the parent halo are rare I Torm en et alJll99^IChen et alJl2003h . so 
we only include subhaloes with masses m O.OlAf. 

We calculate the weak-lensing convergence using Limber's 
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approximation jLimbedflisi lKaiseJll99l . so that the conver- 
gence power spectrum and bispectrum are given by 



and 



I 



diix) VdA(x) 



(14) 



d-iix) VrfA(x)'dA(x)'dA(x) 



(15) 



In the next section we will also use the reduced bispectrum, defined 
as B/3P^ for equilateral triangles. 

We use the standard weak-lensing weight function (see, e.g. 
iBart elmann & Schneider 2001) 



2 dAixs) 



(16) 



where x is the comoving distance, Xs is the comoving distance of 
the source galaxies, ciA(x) is the comoving angular diameter dis- 
tance and Ho is the present-day value of the Hubble parameter. For 
simplicity, we place all source galaxies at redshift 1 . 

We will calculate the covariance of the convergence field by 
adopting a binning scheme. Let P^" denote the power averaged over 
a bin of width Ai centred on the mode h. Our bins are chosen to 
have width Aj/h = 0.3. In this notation, the covariance c an be 
expressed as iCoorav & Shethl2002tlMeiksin & Whitell999h 



An 
hAi 



n 



(17) 



where /sky is the fraction of the sky surveyed, cr^^ is the per com- 
ponent rms intrinsic ellipticity of the sample galaxies, n is the av- 
erage galaxy number density in the survey. The first term is the 
sample variance; the second the shot noise. The last term, T", is 
the trispectrum contribution, arising from non-Gaussianity in the 
convergence field. 

It is likely that a space-based deep imaging survey will of- 
fer the best prospects for probing the small scales where sub- 
structure effects are important. We adort parameters expected for 
the SNAP satellite lensing survey (e.g. ^assev et al. 2003): n = 
100 gal/arcmin^ /sky = 2 per cent, a^^ = 0.34/\/2 = 0.24. 

We use the covariance of the power spectrum to perform a 
least-squares fitting to a fiducial model P"^ by calculating x'^, 



X 



j2{pr-pnc-,\pr~pD' 



(18) 



and identifying la, 2a and 3a contours as those with Ax^ ~ 
2.30, 6.17 and 11.8, respectively, which is appropriate for the two- 
parameter fits we will use. 

We use a ACDM cosmological model with parameters = 
0.3, = 0.7, Qb = 0.05, h = 7, as = 0.9 and the BBKS 
transfer function iBardeen et alll98^ with the shape parameter of 
[Sugivamatl99l . 



3 RESULTS 

Fig. Q compares halo-model calculations of the three- dimensional 
power spectrum with the non-linear fitting formula of ISmith et alj 
(i2003 ). In this figure, different fractions of the halo mass are placed 
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Figure 1. Halo-model calculations of the dimensionless power spectrum 
with different fractions of the halo mass placed in NFW subhaloes. 



in substructure. Adding substructure can be seen to increase power 
at scales fc > 1 Mpc^^. Substructure models predict more power 
than Smith et al. for 10 < fc < 100 Mpc"\ but the model with 10 
per cent of the halo mass in substructure shows better agreement 
at smaller scales (though it should be noted that Smith et al. has 
been extrapolated beyond fc > 40 ft Mpc~^). The power spectrum 
is fairly insensitive to the exponent of the subhalo mass function, 
/i. Reducing /i from 0.9 to 0.7 can increase the power spectrum by 
about 4 per cent at length scales fc > 1 Mpc^^. 

In Fig. |2| we divide the power spectrum into contributions 
from the various terms in equations <A3t - fA71 . On scales fc < 
200 Mpc^^, the substructure contribution to the power spectrum 
is dominated by the density correlation between the smooth com- 
ponent of the parent halo and a subclump, denoted Pac- The cor- 
relation of the density in a single subclump. Pic, or between two 
different subclumps, P2c, is negligible for parameter values of in- 
terest except for very high wavenumbers (where baryonic effects 
are likely to play a role). One should also note that the power of 
the smooth component only, Pss, still dominates all other contribu- 
tions on all scales considered, except for the largest scales where 
the two-halo term, P2h, models the linear power spectrum. 

In Fig. |3| we leave the integral over subhalo mass uncalcu- 
lated for the Pic and Psc terms, so that one can see which subhalo 
masses provide the most power. The figures suggest that most of the 
extra power due to substructure between 10 ^ fc ^ 100 Mpc~^, 
comes from subhaloes of mass 10"' < m < 10^^ Mq. Sub- 
haloes of this mass are likely contained in parent haloes of mass 
> 10^^-10" Af0. Finally, since the weak-lensing efficiency peaks 
for z ~ 0.3-0.5, the conclusion is that the weak-lensing power 
spectrum can probe intermediate-redshift group and cluster mass 
haloes (provided the signal-to-noise permits its measurement as 
discussed below). This should make weak lensing as a probe of sub- 
structure a good complement to studies of strong lensing of quasars, 
which probe substructure in galactic haloes. 

In Fig.|4|we plot a number of power spectra calculated with 
various models. Henceforth, curves titled 'Bullock et al. (2001)' 
and 'Huffenberger & Seljak (2003)' refer to halo-model cal- 
culations without substructure, using the concentration parame- 
ters found in those references. It should be noted that Huffen- 
berger & Seljak obtained their concentration parameter by fitting 
to .Smith etaL L2003) only up fc = 40/i Mpc^V Both Huffen- 
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Figure 2. The dimensionless power spectram separated into the contribu- 
tions from each of the terms in equations 1A3I - IA71 . The order of the Unes 
in the legend matches the order at A; = 10^ Mpc~^. 



Figure 4. Comparison of dark-matter power spectra predicted by various 
models, including substructure, a log-normal distribution of halo concen- 
tration parameters as in equation |9|, and the non-linear fitting formula 
of Smith et al. The order of the lines in the legend matches the order at 
k = 10^ Mpc-i. 
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Figure 3. Contributions from the smooth-subclump term to the power spec- 
trum per logarithmic mass bin. Integrating these curves against d log m 
gives Pac(fc). 



berger & Seljak and Smith et al. have been extrapolated beyond 
k — 40h Mpc^^, and the apparently poor agreement of their 
model with Smith et al. beyond this wavelength should be attributed 
to this extrapolation. 

It is apparent that adding either substructure or the concen- 
tration parameter distribution affects power only on scales k > 
10 Mpc~^. Substructure provides a 10-60 per cent increase in 
power on the scales 10 ^ k ^ 1000 Mpc^^, with the largest 
factor corresponding to the smallest scales. The c-distribution pro- 
vides only a 6-10 per cent increase in power over the same scales. 

It is interesting to not e that while the su bst ructure model 



increa ses power relative to Smith et al. ( 2003) or Bulloc^^l 
|20^|), the concentration parameter of Huffenberser & Selial; 
y003) decreases the relative power. And so, while all models may 
be consistent with Smith et al. at scales ^ 10 Mpc~^, the better 
model should reveal itself as smaller scales are probed in higher- 
resolution simulations and observations. If the halo model is an 
appropriate physical description of mass clustering, then the sub- 



structure model presented here should provide the most robust pre- 
dictions for the power spectrum on small scales. 

In Fig.|5|we show the dimensionless convergence power spec- 
trum predictions for the various models, assuming source galax- 
ies at redshift 1. Again, c-distribution and substructure each pro- 
vide about a 10 per cent increase in power at small scales, though 
the effect of substructure gets larger for th e smallest scales. T he 
halo-model calculations are discrepant with lSmith et alj i2003h at 
the 10-20 per cent level. Most of the contribution to P'^{1) for 
10^ s: / s; 10^ comes from P{k) with 1 fc sC 100 Mpc"\ 
At these scales, the transition between t he one-halo and two-halo 
terms is important. IZehavi et alJ i2003h have proposed a scale- 
dependent correction to the bias parameter and a method to account 
for halo exclusion that may improve the halo model on intermedi- 
ate scales. Takada & Jain ( 2003) employ an approximate correction 
for halo exclusion to obtain a more well-behaved three-point corre- 
lation function on scales ~ 1 Mpc. Given the error bars expected 
from future surveys, it is clear that improved precision in simula- 
tions as well as model predictions is needed. Some additional cos- 
mological models are plotted in Fig.|6| which we will comment 
more on shortly. 

In Figs.0and|8| we plot the predicted three-dimensional bis- 
pectrum and convergence bispectrum. It is interesting to note that, 
though substructure increases the bispectra, it causes a decrease 
in the reduced bispectra. That is, the substructure contribution to 
the power spectrum factor in the denominator exceeds its contribu- 
tion to the numerator in the reduced bispectrum. This is due to the 
fact that the dominant substructure terms are the smooth-subclump 
for the power spectrum and the smooth-smooth-subclump term for 
the bispectrum. The Psc term scales as (1 — /)/, while the _Bssc 
term scales as (1 — /)^/, so the reduced bispectrum scales as 1// 
over most scales. Thus by using both the bispectrum and power 
spectrum, degeneracies between the substructure fraction and other 
cosmological parameters can be broken. We discuss this further in 
Section^ 

Shown in Fig.|9|is the substructure contribution to the bispec- 
trum separated into each of the terms in equations <A10> - 1A15> . On 
the scales of interest, the substructure contributions to the bispec- 
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Figure 5. Vaiious model predictions of dimensionless convergence power 
spectrum. Compares the substructure model and c-distribution to non- 
substructure models. The error bars are the diagonal part of the covariance 
calculated with equation 1171 assuming n = 100 gal/arcmin^, /gjjy = 2 
per cent and ('yf^t)^''^ ~ 0.34. The trispectrum contribution to equation 
1171 is approximated by a halo-model calculation. The order of the lines in 
the legend matches the order at I = 10^. 



2x 10" 




Figure 6. Dimensionless convergence power spectra for models with other 
phenomena that produce effects similar to substructure. The error bars are 
the diagonal part of the covariance calculated with equation 1171 as in Fig.lsl 
The order of the Hnes in the legend matches the order at I = 10^. 

trum are dominated by the smooth- smooth- subclump term, Bssc, 
and the bispectrum for individual subclumps, Bic . The meaning of 
these terms is discussed more completely in the Appendix. It is in- 
teresting that the bispectrum contribution from single subhaloes, 
Bic, provides a larger contribution than the smooth component 
alone, -Bsss, for k ~ 10'' Mpc^^, but on these scales a purely 
gravitational calculation is not likely to be valid. 

Results from the WMAP satellite have constrained the tilt and 
run of the primordial power spectrum. Tilt affects the small scale 
power in a way similar to substructure (cf. Fig.|6}, so we consider 
next the feasibility of weak lensing to separate the effects of sub- 
structure from tilt and run in the primordial power spectrum. The 
results of examining the three parameter {f,nB, dris/d In k) are 
shown in Fig. llOl We use the / = 10 per cent substructure model 



100 r . ■ 

Reduced Bispectrum 

80^ , 

Huffenberger & Seljak (2003) 
"sj ; - Bullock etal. (2001) 
fin - — Substructure Model 

^ i 



■if 




0^ . . ^^^^ 

10 10" 10 

i(l / Mpc) 



Figure 7. Various model predictions of the equilateral bispectrum. The up- 
per panel shows the dimensionless equilateral bispectrum; the lower panel 
is the equilateral reduced bispectrum. 

as the fiducial model, and calculate the for P'^{1) relative to the 
fiducial model using modes in the range 10^ ^ I ^ 10^ The co- 
variance matrix used is that in equation <17t . In the upper panel, 
we assume no running and use Us = 1.00 for the fiducial model. 
In the lower panel, w e fix ria = 0-93 and allow a running spectral 
index, as described in lSoergel et al.'('2003'). to vary about a fiducial 
value of djis/dln k = -0.047 ispergel et aU2003tV Note that we 
do not marginalize over the third parameter in these figures. We 
also apply priors as expected from the upcoming Planck mission: 
fn, = 0.008 and cr^„^/^i^^k = 0.004 lEisenstein. Hu & Tegmar^ 
1998). These figures suggest that with cosmic microwave back- 
ground (CMB) priors, future surveys that can measure P'^ down 
to 0.1-arcmin scales can constrain the substructure mass fraction to 
about ±2 per cent. 



4 DISCUSSION 

We have shown that on small scales, the contribution of CDM sub- 
structure is important in accurate predictions of the power spec- 
trum. In particular, it may help resolve the small scale discrep- 
ancy between TV-body simulation results and the halo model. Our 
model with 10 per cent substructure and a concentration parameter 
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Figure 8. Vaiious model predictions of the equilateral convergence bispec- 
trum. Only the one-halo term is plotted. The upper panel shows the dimen- 
sionless equilateral convergence bispectrum; the lower panel is the equilat- 
eral reduced convergence bispectrum. 
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Figure 9. The one-halo term of the equilateral bispectrum separated into the 
contributions from each of the terms in equations lAlOl - iAlSI . The order 
of the lines in the legend matches the order at A; = 10"^ Mpc~^. 




/ 

Figure 10. Contour plots on two-dimensional planes through -space (no 
marginalization). In the upper panel, the spectral index and substructure 
mass fraction vary with no running spectral index; in the lower panel, the 
spectral index is fixed at ris = 0.93 and a running parameter is varied 
along with the substructure mass fraction. The innermost contour of each 
plot is a one-sigma contour where we have applied prior constraints as ex- 



pected from the Planc k mission! Cj^^ = 0.008 and fj^ 



/d in k 



0.004 



lEisenstein e t alll998h . The other contours are one-, two- and three-sigma 
with no priors. 



distribution fits the small scale power spectrum extrapolated from 
ISmith et'al] i2003h quite well. Higher-resolution A'^-body simula- 
tions are needed to t est our model. 

We agree with Huffenbereer & Seliak (2003) that a concen- 
tration parameter distribution alone canno t resolve the discrepancy 
between the halo model and lSmith et alJ |^003). We have shown, 
however, that the distribution increases power on small scales by 
about 10 per cent, which is comparable to the accuracy of the 
Smith et al. formula. Hence the scatter in the concentration pa- 
rameter needs to be taken into account on scales smaller than 
fc = 10 Mpc-^ 

Our examination of the prospects for parameter con- 
straints are by no means complete. A more careful treat- 
ment must marginalize over a larger parameter space 
X^(f^m, fib, i^i^, o"8, JT-s, dris/d In fc, /, ...), and account for 
survey limitations such as intrinsic alignment of source galaxies, 
the limits imposed by galaxy deblending issues and B-mode 
contamination. Most parameters, like Sim and ag, are not very 
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degenerate with substructure, since they affect the power spectrum 
and bispectrum on all scales, while substructure modifies the 
statistics only on small scales. Such parameters can be well 
constrained from large scale information and probably not confuse 
a measurement of /. Moreover, other observations constrain Sim 
and (78, such as the cosmic microwave background, supemovae 
observations, galaxy clusters and galaxy clustering statistics. 

The situation is similar for neutrino mass and baryon frac- 
tion, though to a somewhat lesser extent. While these parameters 
affect the lensing power spectrum they also produce some sup- 
pression on scales with / < 10* where substructure has little ef- 
fect. This will help us to constrain fi^, 0,^, and / separately. The 
WMAP results suggest that neutrinos are less massive than 0.2 eV 
jSoergel et all EoOS). We have shown that the power spectrum is 
much less sensitive to such neutrinos than substructure at the 10 per 
cent level (Fig.|6|. With regard to baryons, their abundance is well 
constrained from CMB data and big-bang nucleosynthesis. Future 
observations are expected to reveal baryon oscillations in the weak- 
lensing signal which will allow for even tighter constraints on the 
baryon abundance. 

The effects of the parameters cq, fH and a, which denote the 
normalization and slope of the halo concentration and the inner 
slope of the halo profile, can be separated from substructure ef- 
fects with information from the reduced bispectrum. As we noted 
in Sec.|3| substructure produces a proportionately larger increase in 
the power spectrum than in the equilateral bispectrum. The result is 
that substructure amplifies the power spectrum, but attenuates the 
reduced bispectrum. This differs from the behavior of co, a and /3, 
as demonstrated by Takada & Jain ( 2003), who find that these two 
parameters change the power spectrum and reduced bispectrum in 
the same sense. 

Of the parameters considered, the tilt and run of the primor- 
dial power spectrum appear to produce an effect most similar to 
halo substructure. Our analysis in the previous section shows that 
the next generation of weak-lensing surveys will likely provide in- 
teresting Us — f and dns/dlnfc — f constraints. Applying priors 
such as results from the Planck mission to these constraints should 
allow a detection of substructure. We find that a substructure mass 
fraction constraint at the ±2 per cent is possible if observational 
errors can be controlled. 

Numerical simulations indicate that, as a result of tidal dis- 
ruptions in particular, a modified NEW profile for the distribution 
of subhaloes in the parent halo may be more appropriate. As we 
have already pointed out, simulations have not yet provided us with 
a definitive substructure-abundance profile for the centre-most re- 
gions of halos. Even taking t he abund ance to be zero within a core 
radius of order 10 kpc, as in lChen et al. (2003), we would redis- 
tribute only 4 per cent of our subhaloes. So while the magnitude of 
the Pac term may be decreased by this order on 100-Mpc^* scales, 
smaller scales are dominated by the Pic contribution (cf. Eig.|2j, 
which does not depend on the subhalo distribution. 

We thank Gary Bernstein and Ravi Sheth for helpful dis- 
cussions. This work is supported in part by NASA through grant 
NAG5-10924 and NSE through grant AST03-07297. 
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APPENDIX A: HALO-MODEL EXPRESSIONS WITH SUBSTRUCTURE 

We present here th e full expressions for the halo-model power spectrum and bispectrum. These expressions are derived in detail in 
IshethTjai^jlool. 

In the halo-model formalism, the power spectrum is expressed as a sum of one- and two-halo terms: 

=Plh(fc)+P2h(fc). (Al) 

We add substructure by expressing the one-halo term as 

Plh = Pss + Psc + Pic + P2c, (A2) 

where 

.diV(Af) fMA\.^ 



P,,(k) = J AM ^'^^ ' (^-ij U'(k;Ms) (A3) 

represents correlations between the smooth component of a halo. This term looks like the one-halo term of a halo model without substructure, 
but note that the mass of this smooth component is smaller by a fraction /, since Afs = (1 — f)M. Hence, Pss is, for all k, less than the 
one-halo term calculated ignoring substructure. It is up to the remaining substructure terms in equation <A2> to make up for this loss. 

P,,(fc).2 /dM^i^i^C/(fc;Ms)C/.(fc;A/s) /dm^il^§^ (A4) 
J aM p J am p 

denotes correlations between the smooth component and a subclump. The factor of 2 arises because either of the two density elements of the 
two-point function may lie in a subclump. 



dM J dm \ p 

represents correlations where both elements lie in the same subclump and 



P,^^k)^ f dM^^U^ik;M.)\ f dm^^^^^^uik;m 
J dM J dm p 



1 2 

(A6) 



describes correlations between different subclump. Though this expression assumes that subclumps obey a Poisson distribution, we do not 
consider it important since this contribution to the power spectrum is subdominant (see Fig.|2}. 
Finally, 

P..ik) = [ dMi Ml Uik; Ml) I dM, Ml Uik; M.) P(fc; Mi, M.) (A7) 

J dMi p J dM2 p 

is the two-halo term. Note that we have ignored substructure for this contribution; on large scales we trust the approximation of smooth NFW 
haloes with no substructure. 

The one-halo term of the bispectrum has a similar division: 

Plh ~ Bbbb + Esse + Pslc + Bs2c + Pic + P2c + Psc . (A8) 

. dN{M) /Ms^ ^ 



Psss(fci, fca, ka) = J dM \^-y j U{ki;M,) [/(fe; A'/s) U{k3;M,) (A9) 

is the contribution with all three density elements, i.e. all three vertices of the triangle, in the smooth component of a halo. 

2 



Bssc(fci,fc2,fc3) = J U{ki;MB)U(k2;MB)Uc(k3;Ms) 



f dn(m; M) m 

X / dm — V- -—u{k3;m) 

J dm p 

+ cyclic ki permutations (AlO) 
places only one vertex of the triangle in a subclump. 

Bsic(A:i,fc2,fc3) = / dM '^^^y — [/(fci; Afs) Uc{ki- M,) 
J dM p 

f dn{m; M) ( m V 

X / dm ; — u{k2\m)u[k-i\m) 

J dm \p J 

+ cyclic ki permutations (All) 
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has one vertex in the smooth component and the other two in the same subclump. 
J AM p 



X 

1=2 

+ cyclic ki permutations (A 12) 



has one vertex in the smooth component and the other two in different subclumps. 



AN{M) f An{m;M) 



m 



3 



Bic{ki,k2,ki) = J AM ^j^^ J Am — — \~F ) "('^i' ^(^2; m) «(A;3; m) (A13) 
places all three vertices in the same subclump. 

B24k.M,k,) = / AM^^uUk.;M.) f dm /'^^' !^ .(fc,; m,) 

J aM J Anil p 

f An{m2;M) f m2\^ - - 

X / Am2 u(k2;m2)u(ki;m,i) 

J Am2 \ P J 

+ cyclic ki permutations (A14) 
divides the three vertices between two subclumps. 

S3.(fc.,fe,fc3)= / dM^^flU4k.;M.) f dm.^i!i^iMl^.(fc;m.) (A15) 
J J nil p 

places the three vertices in different subclumps. 

This paper has been typeset from a TgX/ LJTjiX file prepared by the author. 
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